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Abstract 



CMS potential for study of the s-channel if ± -boson production via light quark annihilation is in- 
vestigated for large values of tan (3 (— 20 ~ 50) and relatively light charge Higgs boson (Mh = 
200 -7- 400 GeV). An appropriate parameterisation for My (j, ^Y) -distributions of the signal and back- 
ground events is proposed, which provides the determination of the charged Higgs mass and tan (3 
parameter. 
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Introduction 

A charged Higgs boson (if ) appears in many well motivated extensions of the Standard Model (SM) (i.g. in 
the MSSM [[!]]). Two free parameters, M# and tan/3, determine all properties and couplings of the iJ ± -boson. 
The searches for this object were performed in many experiments. LEP collaborations have set lower limits in a 
model independent way on the mass of .H^-boson, M^) > 78.5 GeV for any tan (3 Two experiments at 
the Tevatron, CDF and D0, have performed several searches for H ± They excluded the low (< 1) and high 
(> 40) tan (3 region up to 120 GeV and 160 GeV, respectively. The additional indirect limit can be evaluated from 
low-energy physics, for example from the values of D s — ► tv and B — ► tv decay branching ratios [Q]. 

The discovery potential of i? ± at LHC has been investigated by both ATLAS [||] and CMS [|[ 0] collaborations. 
It was established that for heavy charged Higgs with M(iJ ± ) > m t the most promising channel is top-Higgs 
associated production in two subprocesses || 

gb^tH*, gg^tH^b 

In the present work we consider an additional subprocess of H ± production due to annihilation of the light q<f -pair 
from the initial hadrons (protons): 

qcf — > H , q — d,u,s,c,b 

The i? ± -boson production in this s-channel with subsequent i? ± decay into tb pair was investigated earlier (see 
|To[ pi]]). Here we investigate the CMS potential for study of the s-channel i? ± -boson production via light quark 
annihilation with subsequent H ± decay into a t^v t lepton pair (see also [pT]]): 

qq' -> H ± -> T ± v T 

We perform our calculation for the case of pp-collisions at ^/s = 14 TeV: 

pp -> H±X (1) 

for large values of tan (3 and for four values of charged Higgs mass: 

tan/3 = 50 and M H = 200, 250, 300, 400 GeV (2) 



We use CTEQ5L parameterisation [ |12[ ] of the parton distributions. All estimates of the expected number of events 
are performed for three year low-luminosity run of LHC: 

Cdt = 30 ftT 1 (3) 

For calculation of the signal and background processes we use the event generator TopReX 3.25 ][l3|], while 
the well-known package PYTHIA 6.157 [[llj is explored for modelling of quark and gluon hadronisation. For 
proper simulation of the detector response all generated events are passed through the fast Monte Carlo package 
CMSJET 4.703 _ 



We find that even after application of all appropriate cuts the expected number of the signal events is relatively 
large (~ 10 2 4- 10 3 ) for M H = 200 ~ 400 GeV and tan/3 - 40 4 50. Therefore, this s-channel H ± -boson 
production process allows not only to establish the presence of the signal from -boson, but also makes possible 
to measure the parameters (the mass and tan (3) of the charged Higgs boson. 

1 Signal event generation 

Fig. p]a presents the diagram describing the process under consideration, 

qq 1 — > H , q — d,u^s,c,b (4) 

Note, that Higgs boson couplings to fermions are proportional to the masses of these fermions [[!]]. Therefore, the 
corresponding production cross section has strong dependence on the light quark mass values. In our calculations 
we use so-called "current" values of m q iJT^] : 

m d = 0.009 GeV, m u = 0.005 GeV, m s = 0.150 GeV, 
m c = 1.250 GeV, m b = 4.800 GeV, 
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which are smaller (especially for light d, u, and s quarks) then those values of quark masses used in PYTHIA 
(m d = m u = 300 MeV, etc, see Q). 

We also take into account the radiative QCD corrections to the 2 — > 1 process (Q). In so doing we calculate also 
the NLO processes (2 — > 2) (see diagrams in Fig. ||) 

qq'^H^g, qg^H^q', q' g ~* H ± q (5) 

It is well-known that the consideration of such processes at small k± (where fcj is the transverse momentum of 
the final particle, i? ± , q or g, defined in the centre-of-mass system of the scattering partons) leads to a double- 
counting problem. Indeed, in this region (k± — > 0) the virtual quark, entering the qq' vertex, has very small 
virtuality and may be considered as a on-shell parton. As a result, any 2^2 process (||) can be factorized into 
two subprocesses. The first one is the initial parton (q or g) splitting into two partons 

q -> qg, q' -> q'g, g -> qq- 

The second subprocess is the quark-antiquark annihilation into i? ± , where one quark (antiquark) comes from 
initial hadron, while the second quark appears due to parton splitting. However, such a process was also calculated 
early as the process (Q). 

This problem was considered in details in Jl7|], where the calculations of the complete Q(a s ) QCD corrections 
(reactions (|^)) to the s-channel production process (Q), including QCD-resummation over multiple soft-gluon emis- 
sion was performed. In our numerical calculations we use an approximation, which provides desirable accuracy 
(see jl8| ] for details). It is based on the consideration of the distribution on the charged Higgs transverse momen- 
tum, pj- (H), defined in the initial pp reference frame. In the region of small (H) the basic contribution to Higgs 
production comes from 2 — > 1 process (Q), while the 2 — > 2 process (^) is responsible for Higgs production with 
high p T {H). 

The method of event generation is thus as follows. Firstly, we generate events with H ± -boson production due to 
2 — * 1 process (Q). Any event from this sample will be accepted if the transverse momentum of the charged Higgs 
is smaller than some value po- Then we generate the second sample of events due to 2 — * 2 process ^) with final 
parton transverse momentum fcy > ko. Any event from this second sample will be accepted if pr(H) > po- Thus, 
the resulting (total) sample of i/ ± -boson production events is the sum of two contributions: 

Nipp^H+X) = N ( - 2 ^ 1 ^(pp — > H ± ; pt(H) < po) 

+ N (2 ^{pp^H ± jet;k T >k ,PT(H)>p ) 

We find that the smooth behaviour of the resulting pj- {H )-distribution can be achieved for the following values of 
these parameters: 

fc w 20 GeV, 
and pa = 29.5 GeV, for M H = (200 ~ 400) GeV 

The corresponding pT(#)-distributions are shown in Fig. |[ 

The behaviour of the production cross section for tan (3 — 50 versus Mh is shown in Fig. ^. For tan/3 > 10 
the branching ratio of the i? ± — > t ± j/ t decay is almost independent on tan/3 (see Fig. ||a). Therefore, the tan/3 
dependence of the production cross section for i? ± has a very simple quadratic dependence: 

a{pp H±X) oc tan 2 (6) 

At the same time the Br(iJ ± — » t v) has a strong dependence on the mass of charged Higgs in the region of 
M H = 200 -7- 400 GeV (see Fig. |b) due to opening of the H ± — » tb decay channel. Therefore, a simultane- 
ous measurement of the charged Higgs production cross section and its mass provides a possibility for indirect 
determination of the tan (3 parameter value. 

2 Signal/background separation 

The most important and irreducible background comes from r ± ^ r production via virtual (Drell-Yan) M /± -boson 
exchange (see Fig. [j]b). All other possible sources of background give relatively small contributions and will be 
not considered. 
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It is well-known that due to the different nature of iJ ± // / and W ± f f interactions the final r-leptons, produced 
via H or W boson exchanges have opposite polarisations (see This feature provides an effective way for 
background suppression. In particular, we use the hadronic r — > irv decay mode, where the the two opposite r 
helicity states lead to remarkably different decay pion laboratory momentum distributions. 

For hadronic r decay identification and reconstruction we use algorithm which allows to identify r-jets in the one- 
prong decay mode (see [|(]] for details). This algorithm is based on the fact that the hadronic decays of r-leptons 
from i? ± — > t^v decay are seen as a narrow, low multiplicity "jet" with a large fraction of calorimetric energy 
consisting from a single track. Due to the opposite polarisation of r-leptons produced via HIW boson decays the 
fraction of the total r-jet energy carried away by the charged track relative to the parent r energy is very different. 
For a Vl^-boson mediated decay such a track carries away a significantly smaller fraction of the r energy then for 
the i? ± decay. It can be seen in Fig. |(], which shows the signal (i? ± ) and background (W*) event distribution for 
Rh variable: 

R h = Eih^/EQet) 

The best signal-to-background separation in terms of this variable is achieved for Rq > 0.8. 

To suppress the large W + jet(s) background with W —> tv (r — > hadron + v) decays we also require a strong 
central jet veto: no other jet with Ej- > 20 GeV in \rj\ < 4.5 

To attempt to extract a if ± — * t^v (with r ± — * + X° + v) signal we thus require: 

• one identified r-jet with E T > 50 GeV and \rj\ < 2.0 

• missing transverse energy 'Ej- > 50 GeV 

• jet veto: no other hadronic jets with £q-(j) > 20 GeV in \rj\ < 4.5 

• no other identified objects (leptons, photons) with Et > 10 GeV in \rj\ < 2.4 

• Rh > Ro = 0.8 (to favour iJ ± onto W* as a source of r's) 

The expected number of events for 30 fb _1 and corresponding signal significances (Ns/y/Ns + Nb) after appli- 
cation of all cuts are given in Table [l]. 

3 Parameterisation of M T -distribution 

Since the expected number of the signal and background events is relatively large, Njj ~ (10 2 -j- 10 3 ) and Nb ~ 
10 3 (Table |l|), we can try to determine the parameters of i? ± -boson, namely the mass of charged Higgs (Mh) 
and tan [3. This could be done by fitting the distribution of the transverse mass M T (j, ^t) of the r-jet and the 
missing energy: 

M\ = (E T (j) + ^ T ) 2 - (pt U) + ^t) 2 , (7) 

here E^(j)(p-j(j)) is the transverse energy (momentum) of the r-jet and is the missing transverse energy in 
the event. 

The r-jet identification and reconstruction algorithm does not provide reconstruction of the full r-lepton momen- 
tum due to the undetected neutrino from r ± — > ft, ± (/i°)^ T decays. As a result, the well-known sharp two-body 
decay Jacobian peak in mj- (r, v) -distribution transforms into a wide bump (see ||(], and Fig. |]) in the observable 

M T (j,$ T ). 

The form of this curve results from the "convolution" of the theoretical m-r(r, ^-distribution (where r and v are 
produced in processes (0-^)) and the "fragmentation" (or "decay") of the produced r-lepton into the observable 
hadronic r-jet. This "fragmentation" depends on the experimental device (detector acceptance, resolution, effi- 
ciency, etc) as well as on the r-jet reconstruction algorithm and could not be calculated theoretically. At the same 
time, any appropriate functional form describing this "fragmentation" will provide a suitable parameterisation in 
our case. 

For this "fragmentation" function D T ^j (z) we use a simple parameterisation as follows jl9|]: 

Dr^j(z)ozz a (z -z)\ (8) 

1 ' In what follows the symbol "j" stands for the r-jet 
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where the scaling variable z = Pt{j) Ipt{t) is the ratio of the transverse momentum of the r-jet (pr{j)) to the 
transverse momentum of the parent r-lepton (pt(t))- Contrary to the case of quark fragmentation, the recon- 
structed r-jet momentum may be larger than the momentum of parent r-lepton (see Fig. ^) due to the detector 
resolution and the r-jet reconstruction algorithm. 

We use the D T -^j(z) parameterisation (||) for the fit of the corresponding z-distributions for all considered values 
of charged Higgs mass, namely for to# = 200, 250, 300, and 400 GeV (see Fig. and Table |J). We do not 
need achieve a good fit in the whole region of z. We are interested in the values of z close to unity, because this 
region corresponds to maximal values of MtU, ^t) close to Mfj. Therefore, in what follows we use the set of 
parameters for D T ^j(z) from (|j) given below: 

a = 6.5, A = 3.5, and z = 1.22 (9) 



Then, the distribution of Mt (j, ^Y) could be evaluated by convoluting the m^(r^)-distribution of the r-lepton 
and neutrino (oc 1/ \jMfj — m_i_(rv)) with D T j (z)-fragmentation function (||): 

dM T {3,%T) Jm t /M s Mf / 2 2_ M 2/ M 2 

where Fq is the normalisation factor and Mf is the mass of the charged Higgs boson to be determined from the fit. 

We perform the fit of the Mj(j, ^^-distribution for pure signal events by means of this parameterisation Fs(M-\-, Mf). 
The results of the fit are shown in Fig. [| One can see that the proposed parameterisation (|Io|) provides not only a 
rather good description of the shape of the Mr(J, ^-^-distribution, but also makes possible the determination of 
the fitted parameter Mf, which is very close to the input mass of charged Higgs boson. 

For the Mj (j, ^^-distribution of the background events we use a simple exponential parameterisation as follows: 

F B (M T ) = exp(a + a ± M T + a 2 M T ) (11) 
We get the following values of these parameters (see Fig. |J): 

a = 9.65 ±0.067, a x = -0.0421 ± 0.00072, 

a 2 = 0.000194 ±0.0000261, 6 = 1.769 ± 0.0288 ( } 



4 Signal visibility and measurement of charged Higgs parameters 

In the fitting procedure of the Mj(j, ^-^-distribution of the joint sample of signal and background events we fix 
all parameters in the parameterisations Fs(Mt, Mf) and Fb(Mt) except the corresponding normalisations (Fq 
in Fs and ao in Fb). The mass of the iJ ± -boson (Mf in Fs{Mt, Mf)) is also left as a free parameter to be 
determined by the fit. 

The results of the fit are shown in Fig. [Io| and Table [| One can see that the extracted charged Higgs masses 
(Mf) from the fitting procedure coincide with input values (Mh) within the errors. Therefore, the proposed 
parameterisation of the M-j(j, ^-^-distributions for signal events provides a reasonable way for a determination 
of the charged Higgs mass. 

Then, using the normalisation parameters Fq and ao for the signal and background events we may evaluate the 
corresponding number of events: 

N f s = [ F s (M T ,Mf)dM T , 

JM T , m in 

N f B = [ F B (M T )dM T 

JM T , m in 

These number are given in Table f| One can see that the values of and Ng extracted from the fitting procedure 
are in a good agreement with the expected (generated) numbers Nb and Ns (see Table [l]). 

As a criterion for detection of the signal we use a significance criterion as follows (which corresponds to 99% CL): 

N f 

, 3 >3 (13) 
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Due to our cuts where we require Ej- (J) (St) > 50 GeV we perform the fit of the Mj (j, *Et) -distributions of the 
signal and background events (see Fig. [l0| and Table |]) for Mj values greater than Mj. m i„ = 100 GeV. 

From Table [}]one can see that the significance of the fit for Mji > 300 GeV is too low. To increase the signal-to- 
background ratio we truncate the My distribution, i.e. we repeat the same fitting procedure, but for new value of 



j, which is equal to one half of Mf ] : 



M| 2) - = -M' 1 ' 

T,mzn o J 



(2) 

and determine the new value for M) . This decreases significantly the number of the background events. The 
results of this new fitting procedure are given in Table Q One obtains almost the same results for Mj, but with on 
increased significance (compare Table [3] and Table Q). 

As explained in Section 1, for Mjj > 200 GeV the production cross section for i? ± has the almost quadratic 
dependence on tan/3 (see (^)). Comparing the number of signal events extracted from the fit (Ng) with that 
expected (Ns) in the MSSM scenario, we could determine the tan (3 parameter by means of simple equation: 



^-^ NsiM^i^^y (14) 

where Ns(Mf,ta.n(3 = 50) is the number of expected events with i? ± -boson, generated with Mh = Mj and 
tan (3 = 50 (our default parameters, see (En). The corresponding uncertainty is evaluated as follows: 



S {tan(3f)= t ^l^ S 2 N+dli (15) 

where Sn is the relative error due to F$ parameterisation. A second relative error (5m ) is due to variation of the 
cross section under Mh variation within its own errors: 

6 M = — \<t(Mj - AM) - a(M f + AM)\ 

The values of fitted Mj for several input tan/? values are given in Table |] and Fig. [TJ]. Naturally, decreasing 
tan (3 leads to decreasing the number of signal events. As a result the error in Mj is increased. The corresponding 
extracted values of tan j3 (for input Mh — 200 GeV) are given in the Table || and Fig. 



Using the criterion ( 13 ) we evaluate the area in the {Mh x tan /3)-MSSM plot which could be explored with this 
process of s-channel _ff ± -boson production followed by a decay to r ± ^ T and a hadronic r decays. This region 
(the upper left corner) is shown in Fig. O. 



5 Conclusion 

We investigate the CMS potential for study of H ± -boson production via s-channel annihilation of light quarks. The 
study is made for large value of tan (3{= 20 -^50) and relatively light charge Higgs boson (Mh = 200-^400 GeV) 
where the method is promising. Simple cuts are proposed for signal-to-background separation. After application 
of these cuts a relatively large number of signal events (Ns ~ 10 2 4- 10 3 ) may be expected. Therefore for this 
region of i/ ± -boson parameter space the study of i? ± -boson production is possible with good a significance. 

We find appropriate parameterisations for My(j, ^-7) -distributions of the signal and background events, which 
allow to determine the mass of charged Higgs, which is very close to the input values of Mh- Comparing the 
number of signal events evaluated from the fit, with the expected ones in the framework of the MSSM, we can 
determine tan (3 with a reasonable accuracy. Using a standard significance criterion we determine region in the 
(Mh x tan (3) parameter plot, where this method could be applied. 

Finally, we conclude that the proposed subprocess of s-channel for the iJ ± -boson production does provide a good 
possibility for detecting a relatively light charged Higgs boson and the measurement of its mass and tan j3, 
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Table 1: The total number of the signal and background events after application of all cuts. Signal events are 
generated for tan = 50 and four values of i/ ± -boson mass. The integrated luminosity is Li nt = 30 fb _1 . 



Mh, input 


200 


250 


300 


400 


bkg 


1 v ev 


1627 


344 


129 


35 


1756 


N s 
^/N. q +N B 


28 


7.5 


3.0 


0.83 





Table 2: The parameters resulted from the fit to "fragmentation" function D T ^j(z) from ([$]). The fit was per- 
formed for four i/ ± -boson mass values and in the region of 0.65 < z < 1.2 (z — Pt(J) /.Pt (?"))• 



M H (GeV) 


a 


A 


zq (fixed) 


X z /N 


200 


6.9 ±0.4 


3.35 ±0.17 


1.22 


42. /23 


250 


6.5 ± 1.0 


3.09 ±0.32 


1.22 


3.2/23 


300 


6.6 ± 1.8 


3.49 ±0.61 


1.22 


1.8/23 


400 


6.1 ±4.2 


3.26 ± 1.44 


1.22 


0.23/23 



Table 3: Results of the fit of joint signal and background events Mj(j, ^t) -distribution for M-j > 100 GeV. 



M H 


M f (GeV) 


N B 


iV s 


y/N(,+N* 


200 


202. ±2.1 


1694 


1444 


25.8 


250 


256. ± 9.2 


1694 


250 


5.9 


300 


305. ±20. 


1694 


115 


2.7 


400 


392. ±42. 


1694 


41 


1.0 



Table 4: Results of the same fit as in Table |3|, but for Mt > Mf/2, where Mf is result of previous fit. 



M H 


-iV^T ,min 


Mf (GeV) 


N B 


N f s 




\/ N L+ Nf s 


200 


> 100 


202. ±2.1 


1694 


1444 


25.8 


250 


> 125 


256. ±9.4 


899 


231 


6.9 


300 


> 150 


300. ± 19. 


506 


97 


4.0 


400 


> 200 


391. ±43. 


245 


32 


1.9 



Table 5: Results of the fit of joint signal and background events Mt(j, ^redistribution for Mh = 200 GeV and 
several input values of tan (3. 



tan (3 


50 


40 


30 


20 


15 


Mf 


201 ±2 


203 ±3 


205 ±5 


212 ± 13 


222 ±28 


tan/3/ 


48.3 ±2.6 


39.3 ±4.1 


31.3 ±5.2 


19.8 ±7.6 


16.2 ± 10.6 
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Figure 1 : Diagrams describing light quarks (q = d,u, s, c, b) annihilation into r ± v T pair via charged i? ± -boson (a) 
and -boson (b) exchange. 
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qq' ->■ 
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q — »- 



-q' 



il ± 



Figure 2: Diagrams describing NLO corrections to 2 — ► 1 process of iJ ± -boson production (q = d, u, s, c, b) 
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Figure 3: vr (H )-distribution of iJ ± -boson production in LO subprocess (Q) (solid histogram) and in NLO 
subprocess (gj) (dashed histogram). The vertical line corresponds to parameter value po — 29.5 GeV (see text). 
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Figure 4: Behaviour of the total production cross section for i? ± -boson in pp-collisions (reaction ([!])) at y/s = 
14 TeV and tan/3 = 50 versus charged Higgs boson mass, Mh (solid curve). The dashed curve represents the 
same cross section times branching fraction to iJ ± — > t^u t . 
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-b 




\ / Br(H + ^tb) 




/ \ Br(H + -» t + v) 


", , , , 1 , , 




150 


200 250 300 
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Figure 5: The behaviour of the branching fractions of i7 ± -boson decays into t^v t and tb pairs versus tan/3 (a) 
and Mjj (6). 



no cuts 



jet veto + p^x-jet) > 50 GeV 




Figure 6: The number of events distribution versus Rh = E(h ± )/E(jet. The solid (dashed) histograms corre- 
spond to the background (signal) events. The left histogram is before any cuts, while the right one corresponds to 
requirements of "jet-veto" and _E-r(T-jet), i^y) > 50 GeV. 
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Figure 7: Signal events distribution as a function of the ratio of the transverse momentum of r-jet (pT(i~-jet)) 
to transverse momentum of parent r-lepton (pt(t)). The curves are results of the fit to "fragmentation" function 
D T .,(-.: from (g). 
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Figure 8: Mt(j, ^?t) -distribution for the signal events for four input values of i? ± -boson mass. The curves are 
results of the fit to Fg(Mr,Mf) parameterisation from <JTo|). 
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Figure 9: The result of the fit of My (j, ^Y) -distribution for the background events. 
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Figure 10: The results of the fit of joint signal and background Mt(j, ^^-distribution by the sum of F$ and Fb 
parameterisations. 
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Figure 12: The allowed region in (Mjj x tan j3) plot, which will be available for investigation in the of s-channel 
7J ± -boson production (the left upper corner). 
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